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Summary  of  Research.  Results 


This  description  of  results  from  the  project  covers  research  in  two  areas: 

1.  Wave  propagation  for  equations  describing  elastoplastic  deformations  of  granular  ma¬ 
terials. 

2.  Hyperbolic  conservation  laws. 


1.  PLASTIC  FLOW 

This  part  of  the  project  was  carried  out  in  collaboration  with  David  Schaelfer  of  Duke 
University.  Our  research  focused  on  issues  of  wave  propagation  in  the  equations  of  motion  of 
elastoplasticity.  The  partial  differential  equations  are  derived  from  conservation  of  mass  and 
momentum,  augmented  by  constitutive  laws  that  relate  the  dependent  variables  algebraically. 
The  constitutive  laws  considered  describe  properties  of  granular  materials. 

Scale-invariant  problems. 

We  used  scale-invariant  initial- value  problems  (i)  to  clarify  an  example  of  nonuniqueness 
in  multidimensional  nonassociative  plasticity  [1]  and  (ii)  to  construct  solutions  for  certain 
model  problems  in  plasticity  [3]. 

Shear  bands 

References  [5]  and  [4]  solve  the  Riemann  problem  for  a  system  based  on  antiplane  shear¬ 
ing,  for  cases  in  which  strong  loading  forces  a  shear  band  to  form.  The  shear  band  is  treated 
as  a  discontinuity  in  the  material  subject  to  an  ad  hoc  jump  condition.  Mathematically  the 
solution  with  a  shear  band  is  interesting  because  (i)  the  Riemann  problem  does  not  admit 
a  scale-invariant  solution  and  (ii)  the  generalized  solution  leads  to  a  natural  free  boundary 
problem  for  the  wave  equation.  Reference  [5]  analyzed  the  small-time  asymptotics  of  the 
problem;  [4]  proves  a  rigorous  existence  theorem.  The  latter  was  surprisingly  tricky,  appar¬ 
ently  requiring  the  Nash-Moser  theorem.  These  analytic  results  have  played  a  major  role  in 
guiding  computations  with  shear  bands  undertaken  by  my  colleague  Xavier  Garaizar. 

Incrementally  nonlinear  constitutive  laws 

References  [8]  and  [9]  begin  the  theoretical  analysis  of  wave  propagation  with  rate- 
independent,  incrementally  nonlinear  constitutive  laws,  a  class  which  contains  both  yield- 
vertex  models  and  hypoplastic  models.  Physically  these  are  important  in  improving  the  ac¬ 
curacy  of  modeling,  especially  in  the  presence  of  large  stress  rotations,  and  mathematically 
they  are  interesting  because  the  resulting  partial  differential  equations  are  fully  nonlinear. 
These  papers  demonstrate  that,  for  fully  nonlinear  equations,  there  is  a  subtle  relationship 
between  linear  well-posedness  (i.e.,  hyperbolicity)  and  nonlinear  well-posedness  (existence 
of  a  unique  solution  of  initial  value  problems  and  continuous  dependence  on  the  data).  In 
reference  [8],  we  analyze  the  structure  of  discontinuous  solutions  and  identify  a  class  of  2x2 
model  systems,  for  which  well-posedness  of  the  Riemann  problem  can  be  identified  with  the 


property  of  linear  hyperbolicity.  In  reference  [9]  we  present  examples  to  show  that  such 
an  identification  is  not  possible  for  all  models.  In  this  paper  we  also  study  solutions  of 
the  regularized  equations,  deriving  the  structure  anticipated  in  [8].  Reference  [12]  studies 
numerically  the  blow-up  of  the  solution  in  an  ill-posed  case.  In  recent  work  [11],  we  are 
extending  the  analysis  to  the  full  5x5  system  for  hypoplastic  plane  waves  in  two  dimensions. 
A  preliminary  account  of  this  work  appears  in  [10]. 

2.  HYPERBOLIC  CONSERVATION  LAWS 

There  are  classes  of  nonstrictly  hyperbolic  systems  for  which  the  characteristic  speeds 
are  real,  but  coincide  along  a  curve.  In  [2],  we  give  a  local  analysis  and  classification  of 
such  equations,  and  show  how  new  types  of  shocks,  known  as  singular  shocks,  are  a  crucial 
part  of  solving  Riemann  problems.  We  give  a  detailed  asymptotic  analysis  of  corresponding 
solutions  of  the  regularized,  parabolic,  equations.  These  solutions  blow  up  at  a  single  point 
as  the  dissipation  approaches  zero. 

In  [7],  based  on  the  masters  project  of  my  student,  I  show  how  undercompressive  shocks 
are  the  limit  of  travelling  wave  solutions  of  a  scalar  hyperbolic  conservation  law  with  added 
dissipative  and  dispersive  terms.  In  this  paper,  the  Riemann  problem  is  solved  using  the  new 
shock  waves,  and  numerical  results  are  given  that  show  clearly  the  stability  of  the  travelling 
waves.  The  analysis  depends  on  new  results  concerning  saddle-to-saddle  trajectories  for 
cubic  vector  fields  in  the  plane. 

A  much  more  complicated  analysis  was  performed  earlier  [6]  for  a  p-system  of  mixed 
type,  forming  the  thesis  of  Yadong  Yang.  The  paper  demonstrates  precisely  the  regime  of 
initial  data  for  which  admissible  solutions  are  nonunique. 
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